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5.6 Nonlinear Flow and Conservation Laws 

Nature is nonlinear. The coefficients in the equation depend on the solution u. In 
place of Ut = cux we will study Ut + uu^ = and more generally Ut + f{u)x = 0. 
These are "conservation laws" and the conserved quantity is the integral of u. 

The first part of this book emphasized the balance equation: forces balance and 
currents balance. For steady flow this was Kirchhoff's Current Law: flow in equals 
flow out. The net flow was zero. Now the flow is unsteady — the "mass inside" is 
changing. So a new d/dt term will enter the conservation law. 

There is "flux" through the boundaries. In words, the rate of change of mass 
inside a region equals that incoming flux. For an interval [a,b], the incoming 
flux is the difference in fluxes at the endpoints a and b: 

Integral form — / u{x,t) dx = f{u{a,t)) — f{u{b,t)) . (1) 

at J a 

In applications, u can be a density (of cars along a highway). The integral of u gives 
the mass (number of cars) between a and h. This number changes with time, 
flow in at point a and out at point h. The flux is density u times velocity v. 

The integral form is fundamental. We can get a differential form by allowing h to 
approach a. Suppose b — a = Ax. If u{x,t) is a smooth function, its integral over a 
distance Ax will have leading term Axu{a,t). So if we divide equation (1) by Ax, 
the limit as Ax approaches zero is du/dt = —df{u)/dx: 




When f{u) = density u times velocity v{u), we can solve this single conservation 
law. For traffic flow, the velocity v{u) can be measured (it will decrease as density 
increases). In gas dynamics there are also conservation laws for momentum and 
energy. The velocity v becomes another unknown, along with the pressure p. The 
Euler equations for gas dynamics in one space dimension include two additional 
equations: 

d d 

Conservation of momentum —(uv) + —iuv'^ +p) = (3) 

ot ox 



Conservation of energy TriE) + —(Ev + Ep) = . (4) 

at ox 

Systems of conservation laws are more complicated, but our scalar equation (2) al- 
ready has the possibility of shocks. A shock is a discontinuity in the solution m(x, t), 
where the differential form breaks down and we need the integral form (1). 
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The other outstanding example, together with traffic flow, is Burger's equation, 
for u = velocity. The flux f{u) is 2'^^- The "inviscid" form has no Uxx- 

„ , . du d f du du ^ 

Burger s equation Yt ^ diKj ) = Tt ^ = ^ ' 

When both the density and velocity are unknowns, these examples combine into 
conservation of mass and conservation of momentum. Typically we change density to 
p. For small disturbances of a uniform density po? we could linearize the conservation 
laws and reach the wave equation (Problem ). But the Euler and Navier-Stokes 
equations are truly nonlinear, and we begin the task of solving them. 

We will approach conservation laws (and these examples) in three ways: 

1. By following characteristics until trouble arrives: they separate or collide 

2. By a special formula ( ) 

3. By flnite difference and flnite volume methods, which are the practical choice. 



Characteristics 

The one-way wave equation ut = cux is solved by ^(x, t) = u{x + ct, 0). Every initial 
value Mo is carried along a characteristic line x + ct = xq. Those lines are parallel 
when the velocity c is a constant. 

The conservation law Ut = +uux = will be solved by u{x,t) = u{x — ut,0). 
Every initial value uq = u{xo, 0) is carried along a characteristic line x — uqI = Xq. 
Those lines are not parallel because their slopes depend on the initial value Uq. 

Notice that the formula u{x,t) = u{x — ut,0) involves u on both sides. It gives 
the solution "implicitly." If the initial function is u{x, 0) = 1 — x, for example, the 
formula must be solved for u: 

1 — X 

M = 1 — (x — ut) gives (1 — t)u = 1 — X and u = - — - . (5) 

This does solve Burger's equation, since the time derivative Ut = {I — x)/{l — t)^ is 
equal to —uux. The characteristic lines (with different slopes) can meet. This is an 
extreme example, where all characteristics meet at the same point: 

X — UQt = Xq or X — (1 — Xo)t = Xo which goes through x = 1, t = 1 (6) 

You see how the solution u = (1 — x)/(l — t) becomes 0/0 at that point x = 1, t = 1. 
Beyond their meeting point, the characteristics cannot completely decide u{x,t). 

A more fundamental example is the Riemann problem, which starts from two 
constant values u = A and u = B. Everything depends on whether A > B or A < B. 
On the left side of Figure 5.13, with A > B, the characteristics meet. On the right 
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side, with A < B, the characteristics separate. Both cases present a new (nonhnear) 
problem, when we don't have a single characteristic that is safely carrying the correct 
initial value to the point. This Riemann problem has two characteristics through the 
point, or none: 

Shock Characteristics collide (light goes red: speed drops from 60 to 0) 
Fan Characteristics separate (light goes green: speed up from to 60) 

The problem is how to connect m = 60 to m = 0, when the characteristics don't give 
the answer. A shock will be sharp breaking (drivers only see the car ahead in this 
model). A fan will be gradual acceleration. 

TO DO... 

Figure 5.13: A shock when characteristics collide, a fan when they separate. 

For the conservation law ut + f{u)x = 0, the characteristics are x — f '{uo)t = xq. 
That line has the right slope to carry the constant value u = uq: 

-uixo + St,t) = S — + — = when S = f'iu). (7) 

The solution until trouble arrives is u{x, t) = u{x — f'{u)t, 0). 



Shocks 

After trouble arrives, it will be the integral form that guides the choice of the correct 
solution u. If there is a jump in u [a shock), that integral from tells where the jump 
must occur. Suppose u has different values ul and ur at points xl and xr on the 
left and right sides of the shock: 

d 

Integral form — / udx + /{ur) — /{ul) = . (8) 

dt J^^ 

If the position of the shock is x = X{t), we take xl and xr very close to X. The 
values of u{x,t) inside the integral are close to the constants ul and ur: 

^ [{X - Xl) Ul + (xr - X) Ur] + f{uR) - f{uL) ~ . 
This gives the speed s = dX/dt of the shock curve: 



Jump condition s Ul — s ur + /{ur) — /{ul) = 

U T A fiUR)-f{UL) [/] 

shock speed = = . (9) 

ur-Ul [u] 
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For the Riemann problem, the left and right values ul and m/j will be constants A 
and B. The shock speed s is the ratio between the jump [/] = f{B) — f{A) and 
the jump [u] = B — A. Since this ration gives a constant slope, the shock line is 
straight. For other problems, the characteristics are carrying different values of u into 
the shock. So the shock speed s is not constant and the shock lines is curved. 

The shock gives the solution when characteristics collide. With f{u) = in 
Burger's equation, the shock speed is halfway between ul and ur: 

Burger's equation Shock speed s = = -{ur + ul) ■ (10) 

2 ur-ul 2 

The Riemann problem has ul = A and ur = B, and s is their average. Figure 5.14 
shows how the integral form of Burger's equation is solved by the right placement of 
the shock. 



Fans 

You might expect a similar picture (just flipped) when A < B. Wrong. The integral 
form is still satisfied, but it is also satisfied by a fan. The choice between shock and 
fan is made by the "entropy condition" that as t increases, characteristics must go 
into the shock. The wave speed is faster than the shock speed on the left, and slower 
on the right: 

Entropy condition f'{u) > s > f '{ur)^^^ (11) 

Since Burger's equation has f'{u) = u, it only has shocks when ul is larger than ur. 
In the Riemann problem that means A > B. In the opposite case, the smaller value 
Ul = A has to be connected to ur = B by the fan in Figure 5.14: 

X 

Fan (or rarefaction) ^ ~ ~t At < x < Bt . (12) 

use fig 6.28 p. 592 of lAM (reverse left and right figs) 

Figure 5.14: Characteristics collide in a shock and separate in a fan. 

Notice especially that in the traffic flow problem, the velocity v{u) decreases as 
the density u increases. A good model is linear between v = f max at zero density 
and f = at maximum density. Then the flux f{u) = uv{u) is a downward parabola 
(concave instead of Burger's convex m^/2): 

TYafRc^speed .(.) = ..ax (l - ^) and = -max - ^) • (13) 

Typical values for a single lane of traffic show a maximum flux of / = 1600 vehicles 
per hour, when the density is m = 80 vehicles per mile. This maximum flow rate 
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is attained when the velocity f /u is v = 20 miles per hour! Small comfort at that 
speed, to know that other cars are getting somewhere too. 

Problems and compute the solution when a light goes red (shock trav- 
els backward) and when a light goes green (fan moves forward). Please look at the 
figures, to see how the vehicle trajectories are entirely different form the characteris- 
tics. 

A driver keeps adjusting the density to stay safely behind the car in front. (Hitting 
the car would give u < 0.) We all recognize the frustration of braking and accelerating 
from a series of shocks and fans. This traffic crawl happens when the green light is 
too short for the shock to make it through. 

A Solution Formula for Burger's Equation 

Let me comment on three nonlinear equations. They are useful models, quite special 
because each one has an exact solution formula: 

Conservation law Ut + u 

Burger's with viscosity Ut + uu,j. 
Korteweg-de Vries Ut + uUx 




The conservation law can develop shocks. This won't happen in the second equation 
because the u^x viscosity term prevents it. That term can stay small when the solution 
is smooth, but it dominates when a wave is about to break. The profile is steep but 
it stays smooth. 

As starting function for the conservation law, I will pick a point source: u{x, 0) = 
5{x). We can guess a solution with a shock, and check the jump condition and entropy 
condition. Then we find an exact formula when u u^x is included, by a neat change of 
variables that produces ht = vhxx- When we let — 0, the limiting formula solves 
the conservation law — and we can check that the following solution is correct. 

Solution with u{x,Q) = 5{x) When m(x,0) jumps upward, we expect a fan. 
When it drops we expect a shock. The delta function is an extreme case (very big 
jumps up and down, very close together!). So we look for a shock curve x = X{t) 
immediately in front of a fan! 

X 

Expected solution u{x,t) = — for < x < X{f)\ otherwise m = 0. (14) 

The total mass at the start is j 5{x) dx = 1. This never changes, and already that 
locates the shock position X{t): 

X X'^ ^ 
Mass at time t = \ - dt = — = 1 so X{t) = V2t . (15) 
Jo t 2t 
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Does the drop in m, from X/t = y/2i/t to zero, satisfy the jump condition? 



Shock speed s 



dX 
Itt 



2v^ 



equals 



Jump [mV2] _ XV2t^ 
Jump [u] X/t 



2t 



The entropy condition ul > s > ur = is also satisfied, and the solution ( ) looks 
good. It is good, but because of the delta function we check it another way. 

Begin with Ut + uu.j. = uu^x, and solve that equation exactly. If u{x) is dU/dx, 



then integrating our equation gives Ut + ^U^ = i^Uxx- The initial value Uq{x) is now 
a step function. Then the great change of variables U = — 2i/log/i produces the heat 

equation ht = yhxx (Problem ). The initial value becomes /i(a:, 0) = ^-^o{x)/2u^ 

Section 5.4 found the solution to the heat equation Ut = Uxx from any starting function 
h{x, 0) and we just change t to ut: 



U{x,t) = -2u\og h{x,t) = -2z/log 



y/Anut 



-Uo{y)/2u-{x-y)^/iut 



dy 



(16) 



It doesn't look easy to let z/ ^ 0, but it can be done. That exponential has the 
form e"^'^^'^)/^'^. This is largest when B is smallest. An asymptotic method called 
"steepest descent" shows that as u —>■ 0, the bracketed quantity in (16) approaches 
ce~^~rnin/^'^. Taking its logarithm and multiplying by — 2z/, (16) becomes U = -Brnin 
in the limit: 



lim U (x, t) = -Brnin = min 

u^O y 



Uo{y) + -{x - y)' 



(17) 



This is the solution formula for Ut + ^U^ = 0. Its derivative u = Ux solves the 
conservation law Ut + uux = 0. By including the viscosity uuxx with i/ — 0, we are 
finding the u{x,t) that satisfies the jump condition and the entropy condition. 

Example Starting from u{x, 0) = S{x), its integral Uq is a step function. The minimum 
of B is either at y = x or at ?/ = 0. Check each case: 



U{t,x) = B, 



mm 



min, 



(y<0) {x-yf 

1 (y>0) ^ 2t 



for X < 
xV2t for < X < y2t 

1 for X > V2t 



The result u = dU/dx is or x/t or 0. This agrees with our guess in equation ( ) — a 
fan rising from and a shock back to at x = \/2t. 



